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Abstract 

We study percolation as a critical phenomenon on a multifractal support. The scaling expo- 
nents of the the infinite cluster size (/3 exponent) and the fractal dimension of the percolation 
cluster (df) are quantities that seem do not depend on local anisotropics. These two quantities 
have the same value as in the standard percolation in regular bidimensional lattices. On the 
other side, the scaling of the correlation length (v exponent) unfolds new universality classes 
due to the local anisotropy of the critical percolation cluster. We use two critical exponents 
v according to the percolation criterion for crossing the lattice in either direction or in both 
directions. Moreover v is related to a parameter that characterizes the stretching of the blocks 
forming the tilling of the multifractal. 

PACS numbers: 64.60.Ak, 61.43.Hv, 05.40.-a, 64.60.Fr 
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1 - INTRODUCTION 



One remarkable feature of critical phenomena is the concept of universality, that means, 
the fact that critical exponents are independent of the particular type of lattice and of 
nature of the physical system. Transitions in systems so diverse as gases, magnets or 
metalic alloys, can be characterized by the same critical indices. In some sense this can 
be explained as a consequence of the equivalence of the local simmetry of these different 
systems jl| . Non universal behavior is a rare event observed only in special situations and 
non trivial models. 

In magnetic systems the critical exponents, in general, depend only on the dimension 
of the space and on the dimension of the spins. The variation of critical exponents for 
fixed dimension of the lattice and fixed dimension of the spin is found in non trivial cases 
in which there is some types of competing interactions, frustration or symmetry breaking 



In this work we report results obtained for a complex model which is universal with 
respect to some exponents but that shows non-universality with respect to the correlation 
length exponent v. We deal with percolation phenomena in a multifractal support im- 
mersed in a 2D space. It has been shown that percolation corresponds to the q = 1 state 
of Potts model [4]. In this perspective we expect universal behavior for percolation in 2 
dimensions, for any kind of lattice. The multifractal support we have proposed has a new 
local lattice structure: the local connectivity and the local anisotropy vary from region 
to region. The strength of the anisotropy change of the multifractal support is related 
to a single parameter p. The exponent v of the critical percolation cluster varies with p 
showing an accentuated sensibility to the local anisotropy of the support structure. On 
the other side, the percolation cluster is a fractal with a fractal dimension, df, that is 
identical to the one of the standard percolation cluster. The same occurs with the critical 
exponent f3. That means, if we take into consideration only df and j3 we could say that 
the universality class holds for this problem. 

A motivation for this work also comes from the problem of direct percolation in which 
the anisotropy of the model is responsible for anisotropic clusters and for different correla- 
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tion lengths in parallel and perpendicular directions. In this work we observe the build up 
of percolation clusters that are very elongated for small values of p, showing anisotropic 
structures. Do these local anisotropic properties affect the critical exponents? 

In the reference j^, i| an object, Q m f, is developed to study percolation in multifractal 
systems. This object is a natural generalization of the square lattice in which it is possible 
to estimate analytically the full spectrum of fractal dimensions. Besides it exhibits a rich 
distribution of area among the blocks and a non-trivial topology, the number of neighbors 
varies along the object. Moreover its geometric qualities and algorithmic simplicity, it is 
straightforward to study the percolating properties of Q m f- 

In a previous work 5| we have introduced the multifractal, analytically determined its 
spectrum of fractal dimensions, analyzed some properties of the histogram of percolating 
lattices versus occupation probability, and estimated the fractal dimension of the perco- 
lation cluster for several parameters of Q m f- In a further work p we have performed a 
comprehensive study of the determination of the percolation threshold, p c . In particular 
we have explored the relation between p c and the topologic properties of Q m f. In this 
work we deal with the problem of the class of universality of Q m f and its relation with 
symmetry characterization. 

The multifractal object we develop, Q m f, is an intuitive generalization of the square 
lattice. Suppose that in the construction of the square lattice we use the following algo- 
rithm: take a square of size L and cut it symmetrically with vertical and horizontal lines. 
Repeat this process n-times; at the n th step we have a regular lattice with 2 n x 2 n cells. 
The setup algorithm of Q m f is quite similar, the main difference is that we do not cut the 
square in a symmetric way. In the next section we explain in detail this algorithm. 



2 - THE MULTIFRACTAL SUPPORT 

In this section we present the multifractal object Q m f in a double perspective. Initially 
we introduce the object Q m f using its growing algorithm. In a second moment we show 
how percolation properties are studied in such a non-trivial system. 
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2-1 The generating algorithm of Q m f 

We start with a square of linear size L and a parameter < p — - < 1, for integers 
r and s. The first step, n — 1, consists of two sections of the square: a vertical and an 
horizontal. Initially the square is cut in two pieces of area -f- and -7- by a vertical line 
(we use L 2 area units). This process is shown in figure C] (a), where we use p = ^ — |. The 
horizontal cut is shown in figure ^(b). Note that we use the same p. The first partition 
of the square generates four rectangular blocks: the largest one of area (tt^;) 2 , two of area 
, t^ni and the smallest one of area (^-) 2 . 

(r+s) z v r+s' 

The second step, n = 2, is shown in figure E(c) and (d). The same process of vertical 
and horizontal sections is repeated as in step 1 inside each block. In this way Q m f is 
a self-affine object. As observed in the figure at end of step n = 2 there are 2 4 blocks. 
Generically we have after the n* h -step 2 2n blocks. The partition process produces a set of 
blocks with a variety of areas. We call a set of all elements with the same area as a k-set. 
In the case n = 2 there are five fc-sets. 

At the n^-step of the algorithm, the partition of the area A of the square in 1? units 
follows the binomial rule: 

^Zq \s + r/ \s + rj \r + sj 

The number of elements of a fc-set is . 

n 

In reference p we see that as n — > oo each k-set determines a monofractal whose 
dimension is dt = lim t i^ 00 iogC '" s r( « — -. In the same limit n — * oo the ensemble of all 

log (s+r) ^ 

A;-sets engenders the actual multifractal Q m f. 

In figure |2 (a) we show Q m f , as in figured] (p = |), for n = 8. We use the following 
code color: blocks of equal area (same k-set) share the same color. Figure 121(b) is a zoom 
of the inner square part of EI ( a )- We can observe in this picture that it shows a slight 
anisotropy compared with the square lattice. The anisotropy in the figure is not very 
accentuated because p = | is close to the square lattice case, p = 1. 

In the construction process of Q m f we see that as r >> s (or p — > 0) the blocks became 
more stretched in one direction than other. In this way p is a measure of the stretching 



4 



of the blocks. This property will reflect in the anisotropy of the percolation cluster as we 
shall see in the simulations. 



2-2 The percolation in Q m f 

The main subject of this work is the study of the percolation properties of Q m f- To 
perform such a task we develop a percolation algorithm. The percolation algorithm of 
Q m f starts mapping this object into the square lattice. The square lattice should be large 
enough that each line segment of Q m f coincides with a line of the lattice, this condition 
imposes that p is a rational number. In this way all blocks of the multifractal are composed 
by a finite number of cells of the square lattice. To explain the percolation algorithm we 
suppose that Q m f construction is at step n. We proceed the percolation algorithm by 
choosing at random one among the 2 2n blocks of Q m f independent of its size. Once a 
block is chosen all the cells in the square lattice corresponding to this block are considered 
occupied. Each time a block of Q m f is chosen the algorithm check if the occupied cells at 
the underlying square lattice are connected in such a way to form an infinite percolation 
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cluster. The algorithm to check percolation is similar to the one used in 

Figure |3] shows the percolating cluster of Q m f for p = | and a particular realization. 
Others p show a similar aspect, but the blocks became more stretched as p — > 0. In this 
limit the multifractal turns more anisotropic. If we compare this percolating cluster with 
the percolating cluster of square lattice p| we observe a strong anisotropy in the present 
case. Such anisotropy is observed also in the scaling of the correlation length, £, as we 
see in the next section. 

Following the literature |7| we call p the density of occupation of a lattice. Rl is 
the probability that for a site occupation p there exists a contiguous cluster of occupied 
sites which cross completely the square lattice of size L. _p c is the density of occupation 
at the percolation threshold. There are several ways |8| to define Rl, we use two of 
them: R e L is the probability that there exits a cluster crossing either the horizontal or 
the vertical direction, and R b L is the probability that there exits a cluster crossing around 
both directions. At the limit of infinite lattice size R e L and R b L converge to a common 
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value in the case of the square lattice. 



3 - NUMERICAL RESULTS 

In this section we deal with the estimation of critical exponents and fractal dimensions 
of Q m f- Before the numerics we warn to a specific point in this process. Suppose the 
multifractal construction is at step n, the underlying square lattice has the number of 
cells 2" x 2", and the number of blocks is 2 2n , values that are independent of p. The size 
of the underlying lattice, however, is L = (r + s) n , a function of p. Because of this point 
we have to take care in the effect of finite size properties of the multifractal. 

n 

A comprehensive study of the percolation threshold, p c , of Q m f is performed in p, we 
resume this point in the following. The way we find the best estimative of p c is to use, 
for several L, the average value of R e L and R b L . The results for specific values of p are 
summarized in Table /. The results point to a value of p c which slightly increases with p 
and present a strong discontinuity at p = 1 (the square lattice). 

3-1 Fractal dimension of the percolating cluster 

The estimation of the fractal dimension of the percolating cluster, df, is performed 
using the definition: 

ln(M) 

d, = hm^—. (2) 

Where M is the mass of the percolating cluster and L the lattice size. We use in our 
simulations values of L that correspond to 5 < n < 10. The results for several p are 
shown in Table I. The exact value of df in two dimensions is df = || = 1.89583. The 
values obtained are close to 2% of the exact value and are not correlated to p. However, 
this fact is not enough to conclude that Q m f is in the same class of universality of two 
dimensional standard percolation. In the next paragraphs we shall return to this point. 
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3-2 The critical exponent (3 



Percolation in bidimensional spaces shows critical phenomenon close to the critical 
point p c . The critical exponent (3 is defined from the relation: 

R L ~(pc(L)-p c ) p , (3) 

where p c is the exact occupation probability value in contrast to p c (L) which is its finite 
size value. The power- law (J3J) is verified for p c (L) obtained either from R e L or R b L . In the 
first case the formation of a spanning cluster and the power-law are at p > p c . Otherwise, 
if we use R b L , we have to search for a power-law in the vicinity of the critical value for 
p < p c . We use in our estimation Rl = R e L which shows smaller fluctuations than R b L . As 
the numerical estimation of (3 is based in equation @, Rl is a key element of the analysis. 
For Q m f the probability Rl is not a well behaved function of p for low L Actually, 
Rl can show, depending on p, an inflection point at p c in this regime. In the case where 
L — ► oo the scaling of {p c {L) — p c ) recovers the power-law behavior. In this regime we 
find the same (3 of two dimension standard spaces, (3 = ^ = 0.13888. We check in our 
simulations that for n = 9, (3 is around 5% of the exact value of standard percolation. A 
set of values of j3 are shown in Table /. 

Table / summarize the estimated p c , df, and (3 for several p. The values of df, and /3 
fluctuate around the theoretical value of the standard percolation and do not show any 
correlation with p. The main conclusion we take from the data is that df and (3 do not 
depend on p. 

Table I 



s,r) 


(1,1) 


(4,3) 


(3,2) 


(2,1) 


(5,2) 


(3,1) 


(4,1) 


(5,1) 


Pc 


0.5929 


0.5253 


0.5267 


0.5262 


0.5260 


0.5261 


0.5254 


0.5253 


df 


1.895 


1.892 


1.890 


1.900 


1.891 


1.911 


1.902 


1.929 




0.127 


0.135 


0.141 


0.128 


0.131 


0.140 


0.141 


0.118 
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3-3 The critical exponent v 



An important quantity in percolation theory is the correlation length, £. It diverges 
at the critical point p = p c with a critical exponent v and it is related to the size of 
the percolation cluster. In a finite size system, near the critical point, we use this fact 
to approximate £ by the system length L. The exponent v is defined from the scaling 
equation: 



In order to estimate numerically v we plot In {p c {L) — p c ) against ln(L). The slope of this 



from R e L . It is indicated in the figure the linear regression equation of the data, v can also 
give some hint about the symmetry of the percolation cluster. In other words, v informs 
about the isotropy of the percolation cluster, or the dependence of £ with some direction. 

We emphasize that Q m f is not an isotropic object and as a consequence the percolating 
cluster is also anisotropic. In this way there is no reason for the percolating cluster to 
be associated to an unique scaling relationship as in equation (@J). We note that R e L is 
related to the existence of a spanning cluster in either direction, and R L in both directions. 
Actually R b L measures an average of the crossing probability over both directions. This 
average over both directions erase the anisotropic effect of the spanning cluster on the 
statistics. Otherwise, the probability R e L (crossing at either directions) accentuate the 
anisotropy of the percolating cluster. Therefore, we do not expect a same result for v 
using p c (L) obtained from R\ or R\. 

Figure 03 shows v versus p for the two cases: crossing in either directions (diamonds), 
and crossing in both directions (triangles). We call v e and v h these two exponents. The 
values of p are indicated in the figure. Based in this figure we initially conclude: diversly 
from df and (3, the exponent v depends on p. For the case p = 1 (the square lattice) we 
have v e = 1.33 and v h = 1.39, considering the sensitivity of the method, these values are 
in good agreement to the exact value of v — ~ = 1.333. 

An anisotropic percolating cluster has a larger spanning probability in either direction, 
R e L , than in both directions, R\. Therefore, as the anisotropy of the cluster increases, the 




(4) 



curve is v. Figure 0] illustrates this process for p 



| and the case where p c is obtained 
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difference between the two estimations becomes more accentuated. This result is observed 
in figure 03 The same result can be verified using Ais = \v e — v h \ which increases as p — > 0. 
In Table II we show the data of figure El and Aza 

Table II 



s,r) 


(1,1) 


(4,3) 


(3,2) 


(2,1) 


(5,2) 


(3,1) 


(4,1) 


(5,1) 


v e 


1.33 


3.80 


3.25 


2.49 


4.78 


3.48 


4.92 


6.13 


v b 


1.39 


3.84 


3.16 


2.03 


3.53 


2.65 


2.98 


3.41 


Ap 


0.06 


0.04 


0.09 


0.46 


1.25 


0.83 


1.94 


2.72 



It is also observed in figure El that generically v decreases with p. This phenomenon is 
related with the stretching effect that increases as p — > 0. We have discussed this point 
in section §2 in connection with the stretching of the blocks that forms the tilling of the 
mult if r act al. 



4 - FINAL REMARKS 



To summarize we study critical exponents and fractal dimensions of the incipient infi- 
nite cluster constructed on a self-affine multifractal Q m f- The quantities (3 and df do not 
depend on the local symmetries of the percolating cluster. As a consequence these two 
numbers are equal to the ones of the standard (isotropic) two dimensional percolation. 
On the other side, the exponent u, which is related to the correlation length, depends 
on the local symmetry properties of the percolation cluster. There are some points to 
consider in the symmetry properties of Q m f and its isotropy. The first is related with the 
stretching of the blocks of Q m f which increases as p — > 0. Moreover there is the topology 
of Qmf, specially its coordination number, that changes along the object. 

We have in mind in the construction of the multifractal object the modeling of fluid flow 



reservoir heterogeneities 
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geophysical situations [15 


, lie. 


17| 



3£, complex geological structures [lj] and to study other 
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The anisotropy of the percolation cluster is revealed when we compare the percolation 
threshold and v in either direction and in both directions. In fact, R\ performs an average 
over both directions erasing the anisotropic effect of the percolation cluster. Moreover we 
estimate two v exponents, v e and v b , which are associated with R e L and R\. Typically 
the correlation length and the exponents v are larger for the percolation over Q m f than 
in usual percolation. This fact is related to the very stretched blocks in the multifractal 
that can make long range connections in the object. We see that the v exponent is able to 
probe the local anisotropy of the percolation cluster on a multifractal support. The value 
of v that is identical to the standard percolation (p = 1, the symmetric square lattice) 
unfolds into different values indicating different degrees of local anisotropics. 

To conclude, we estimate some quantities for percolation on the multifractal object Q m f 
immersed in two dimensions. The exponent (3 and the fractal dimension of the percolation 
cluster are the same of standard percolation in two dimensions. The exponent v that is 
related with the correlation length is influenced by the anisotropy of Q m f- Depending 
on the parameters that characterize the stretching of the tilling of Q m f, p, several values 
of v appear. In this way the universality class of percolation on Q m f depends also on 
p, the anisotropy of the tilling of the multifractal. This surprising result reveals that in 
this problem the local symmetry of the support governed by p affects strongly the critical 
percolation behavior, leading to non-universality. 
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[1] Shang-Keng Ma, Modern Theory of Critical Phenomena , W. A. Benjamin, Reading, 1976 
[2] Domb, C. and M. S. Green, editors Phase Transitions and Critical Phenomena Academic 
Press London, 1976 



10 



[3] Stanley, H. E. Introduction to Phase Transitions and Critical Phenomena Oxford University 
Press, London, 1971 



[4] 

[5] 

[6] 
[7] 

[8 

[9 
[10 

[11 
[12 

[13 

[14 

[15 

[16 
[17 



F. Wu, Rev. Mod. Phys. 54 235 (1982). 

G. Corso, J. E. Freitas, L. S. Lucena, R. F. Soares, cond- mat /021 2530 , accepted in Phy. 
Rev. E (2003). 



L. S. Lucena, J. E. Freitas, G. Corso, and R. F. Soares, cond-mat/0308291 , accepted in 
Brazilian Journal of Physics (2003). 

D. Stauffer and A. Aharony, Introduction to Percolation Theory, (Taylor & Francis, London, 
1994). 

M. E. J. Newman and R. M. Ziff, Phys. Rev. Lett. 85 4104 (2000). 

L. S. Lucena, L. R. da Silva and S. Roux, Physica A 266 86 (1999). 

J. E. de Freitas, L. S. Lucena and S. Roux, Physica A 266 81 (1999). 

J. E. de Freitas and L. S. Lucena, International Journal of Modern Physics C, 11 1581 

(2000). 

F. Herrmann, A Scaling Medium Representation, a Discussion on Well-logs, Fractals and 
Waves. PhD Thesis, Delft University of Technology, (1997). 

J. Muller, Annales Geophysicae - Atmospheres Hydrospheres and Space Sciences 11 (6): 
525 (1993). 

P. N. Khue, O. Fluseby, A. Saucier, J. Muller, Journal of Physics - Condensed Matter 14 
(9): 2347 (2002). 

R. H. Riedi, Multifractals and wavelets: a potential tool in Geophysics, Proceedings of the 
68th SEG Meeting, New Orleans, Louisiana, USA, (1998). 
P. Hubert, Hydrological Sciences Journal 46 897 (2001). 
S. Sachs, S. Lovejoy, D. Schertzer, Fractals 10 253 (2002). 



FIGURE LEGENDS 
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FIG. 1: The two initial steps in the formation of Q m f- In (a) a vertical line cut the square in 
two pieces according to the ratio p = | . Two horizontal lines cutting the rectangles by the same 
ratio are depicted in (b). The second step is show in (c) and (d), the same process of step 1 is 
repeated for each former block. 



FIG. 2: The object Q m f for the same p = | used in figure ^ but the number of steps is evolved 
until n = 8. Figure (b) shows a zoom of the internal square of figure (a). 



FIG. 3: A view of the percolating cluster for one typical realization of p = A 



3" 



FIG. 4: A graphic of ln(p c (L) — p c ) against InL for p = |. It is indicated in the figure the 
equation of linear regression of the data. The slope of the curve is the exponent v. 



FIG. 5: A graphic of values of the exponents v e (square) and v h (circle) versus the anisotropy 
parameter p. 
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